ON AN ALGEBRO-GEOMETRIC DISCRETIZATION OF 

KP HIERARCHY 
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Abstract. This paper studies a certain completely integrable dis- 
cretization of the KP hierarchy. This was constructed by Gieseker 
in []Giel[ , from certain algebro-geomctric data. This paper has the 
dual aim of showing that this construction is generically invert- 
ible, and obtaining explicit expressions for the flow equations. A 
subsequent article will discuss the Hamiltonian structure of this 
system. 



1. The Relation Between the Discretization and 
Algebro-Geometric Data 

1.1. Description of the system. We assume that and M are 
positive integers such that gcd(A^, M) = 1. 

Gieseker describes a family of lattice equations parameterized by 



{N,M) in ||Giel|| . These are two dimensional generalizations of the 
periodic Toda lattice equations and reduce to the latter when M = 1 
. The solutions of the KP equation coming from Riemann G functions 
can be approximated by solutions of the lattice equations, taking M = 
+ 1 and N ^ oo. 
Here, we are going to define the system in a roundabout way. Con- 
sider the following problem: We are looking for functions \E'(n, m,t) 
where {n,m,t) G Z/A^Z x Z/MZ x C that are almost periodic in the 
two space directions of the lattice, i.e.: 



"^(n + N,m) = a'^in, m) 

(1.1) 

where a, j3 are independent of (n, m). Moreover, we require that "^{n, m+ 
1) (suppressing the time variable t) is expressible in terms of three of 
the \E'(fc,m), more specifically that \E'(?T,,m + 1) is of the form: 
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(1.2) 

*(n, m + 1) = *(n + 1, m) - A{n, m)'i>{n, m) - B{n, m)'^{n - 1, m) 

where A{n,m) and B{n,m) are periodic in both space entries, with 
respective periods N and M. Given such a set of A{n,'m), B{n,'m) 
the presence of a nontrivial solution for \1/ forces an algebraic relation 
between a and /?. There exists a matrix W such that the conditions 
above translate as e ker(iy). To get this W, order ^(n, m) keeping 
the second index more significant than the first (i.e. use the order 
(^(1, 1), ^(2, 1), ... ^(iV, 1); ^(1, 2), ...)). W then becomes an NM by 
NM matrix. We present it in block form with N hy N blocks: 



;i.3) w 



xii) 

On 



On 

-In{2) On 
X{2) -In{S) 



J AT 



Jn 



X{M) 
On 



On On X{M-1) -In{M)_ 

W is in block circulant form. It has two nonzero circulants. Block 
(1, 1) of W is — /3 * /iv(l), and for i ^ I, block is — /Ar(i). Block 
{i + of W is X{i) for all i. (Here, regard i in Z/MZ.) /Ar(i) and 
On represent the N hy N identity and zero matrices respectively. The 
sole purpose of indexing Jjv's is making references possible. X{m) is: 



(1.4) 



X{m) = 



-A{l,m) 
-5(2, m) 




a 



1 

-A{2,m) 
-5(3, m) 









1 
-A{3,m) 1 



-B{l,m)/a 




-B{N,m) -A{N,m) 



X{m) is a circulant matrix, this time with three nonzero circulants. 
The (i,i) entry is —A{i,m), the — 1) entry is —B{i,m) , and the 
(i,i + 1) entry is 1 for all i (Here regard i in Z/A^Z). 

We label entries of with two pairs of numbers. Entry ((n,m), {k,l) ), 
where 1 < n, k < N and 1 < m, / < M, means the entry (n, k) of block 
(m, /) of VF. For instance W{{n, m), {n, m—1)) is —A{n, m—1), whereas 
W{{n, m), (n — 1, m — 1)) is —B{n, m — 1). 
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In order for W"^ = and \I/ be nontrivial, det(iy) should be 0. 
Given a set of A{n, m), B{n, m), this is the defining equation of a plane 
algebraic curve in the variables a and f3. We look at a generic element 
of this family of curves. Such a curve has a certain definite behaviour 
at the loci (3 = 0, P = oo,a = and a = oo which we will see in a 
moment. Normalize such that \E'(0, 0) = 1. Then, apart from finitely 
many points of the curve det(iy) = 0, \E'(n, m) is a rational function 
in a and f3. We realize each m) as sections of certain line bundles 
on the normalization X of the curve. Under this connection, the line 
bundle corresponding to \E'(0, 0) turns out to be of degree g =genus of 
the curve. We view this line bundle, C, as a point of the lifted Jacobian 
parameterizing line bundles of degree g on X. 

The interesting thing about this correspondence is that it is possible 
to move the A's and S's in certain ways keeping the curve fixed. This 
can be seen simply by counting dimensions on either side. On the other 
hand, under the correspondence, these degrees of freedom precisely 
correspond to moving C in the Jacobian. From a given curve and a 
line bundle one can retrieve A's and -B's. 

The discrete KP hierarchy is the set of flows corresponding to moving 
C in linear directions on the Jacobian. Since the curve is fixed for 
these flows, one can immediately deduce that there are many conserved 
quantities: the coefficients of the curve equation. Since the flows on 
the Jacobian are linear, the discrete KP flows have to commute. This 
is a completely integrable system. For a discussion of these matters, 
see [|Giel| . 

1.2. The correspondence. We now look at the implications of det(Ty) 
0. We are going to derive most of our results from information about 
monomials fij{A, B)a^P^ that appear in det{W). First of all, we state 
the correspondence between the algebro-geometric data and the dis- 
crete KP data. X denotes the normalization of the curve det W = 0. 

Theorem 1.1. There is a natural correspondence between the follow- 
ing sets of data: 

1) A generic smooth curve X of genus g which possesses points P, Q 
such that N{P — Q) = div(a), an additional list of points Ri,Si, i = 
1,...,M so that M{P + Q) - J^i^i + Si) = div(/3), where a, (3 are 
meromorphic functions on X ; and a line bundle C of degree g on X 
such that 



;i.5) H^{X, C{{n + m - 1)P + (m - n)Q - ^(i?^ + S,)) 

1=1 
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for all {n, m) . 

2) Generic functions A{n,m), B{n,m) , periodic in the two space di- 
rections with periods N and M . 

The relation between g, N and M is g = {N — 1)M. The'^{n,m) ob- 
tained from = 0, if properly normalized, are holomorphic sections 
of the line bundles C{{n + m)P + (m — n)Q — + ^i))- 

The construction of (2) from (1), and hence the introduction of this 
correspondence is the subject of [|Giel|| . It will be described here later. 



^1 studies the reverse construction as well, from the point of view of 



deformation theory. Our approach is more combinatorial, more in the 
spirit of ||vM-M|] . 

The equation of the plane curve is det(iy) = 0. We want to analyze 
the behaviour of this curve when f3 or a tend to oo. We make some 
linguistic conventions: A "monomial" will mean a full multiplicative 
expression involving a, P, A, B. We distinguish one specific piece of a 
given monomial: the "coefficient" of the monomial is the piece formed 
by A's and -B's, in accordance with our treatment of a, (3 as variables. 

Expand the determinant using all permutations of n letters. A mono- 
mial is said to "appear" in the expansion of det(Vr) if there exists a 
permutation vr of NM letters so that the product associated to tt in 
the expansion of det(iy) is a nonzero multiple of this monomial. 

Lemma 1.1. A monomial appearing in the expansion ofdeiiW) with a 
nonconstant coefficient cannot cancel another monomial with the same 
properties in the evaluation o/det(H^). 

Proof: It will be sufficient to prove that a monomial with these prop- 
erties cannot be associated to more than one permutation. Suppose 
that TT is one permutation that such a monomial m is associated to. 
Remove the rows and columns of W that the A's and -B's in m are on. 
We obtain a certain minor H of W . Assign all other A^B^s the value 
0, and call the new matrix Hq. Clearly, Hq has at most two nonzero 
elements in each row or column. Suppose that there is more than one 
permutation on Hq that picks no zero entries. Then, there exists a 
row where these two permutations differ, hence collectively they pick 
both nonzero entries of that row. It follows that they also have to 
pick different entries on the columns of these entries, and so on. The 
complete set of entries picked by one permutation but not the other 
form at least one closed loop that can be traversed by changing either 
the row or column but not both at one move. But this is impossi- 
ble. Indeed, without loss of generality suppose that the diagonal entry 
((6, a), (6, a)) belongs in such a loop. Then, the only other nonzero 
entry of W in this column, ((6 — l,a + 1), a)) should also belong 
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in the loop. Going one step further, the only other nonzero entry in 
this new row, {{b — 1, a + 1), (6 — 1, a + 1)) should belong in the loop. 
Continued, this list contains ((6 — k.a + k), (b — k,a + k)) for all k, 
which exhausts the diagonal of W before coming back to the initial 
point since gcd(iV, M) = 1. Therefore H — W, which is contradictory 
to the assumption that tt picks at least one A or B. □ 

The lemma says that the list of A,B^s in m determines the associated 
permutation tt uniquely, if such a permutation exists. 

Definition 1.1. We assign degrees d to multiplicative expressions in 
a, /9, A, B as follows: 



{i)d{a) 


= N 


{ii)d{(3) 


= M 


iii)d{A) 


= 1 


{iv)d{B) 


= 2 


{v)d{c) - 


= 0,c e C 



and the degree of a product is the sum of the degrees. 

The following lemma suggests that this degree assignment is natural: 

Lemma 1.2. If m is a nonzero monomial appearing in det{W), then 
d{m) = NM. 

Proof: According to the definition, d{{n, 1), (n, 1)) = M, and d{{n, m), (n, m)) = 
for m 7^ 1. Observe that, for nonzero entries of an X{i): 

(1.7) d{{ni, m), (n2, m - 1)) = ni - n2 + 1 

Suppose that m is associated to the permutation tt. Also, suppose 
that the number of /3's that tt picks is k. Let be the multiset (i.e., a 
collection of elements, where elements may be listed more than once) of 
77-1 such that vr((ni, m)) = (n2, m — 1) for some value of m and for some 
77-2. Let 5*2 be the multiset of n2's that appear on the right hand side 
of such an equation for some Ui. If there are t copies of Ui in 5*1, then 
7r(ni, rn) = (n2, m — 1) for t values of m. But for the remaining N — t 
values of m, the only possibility that remains is 7r(ni,m) = (ni,m). 
Now, this implies that ni appears on the right hand sides oi N — t 
equations of the latter form, but then ui satisfies an equation of the 
form Tr{{k, m + 1)) = (ni, m) for precisely the remaining t values of m. 
This implies: 
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(1.8) ^1 = ^2 

Using this and the previous degree calculations: 

d{m) = ^ ni- ^ n2 + M{N ~ k) + Mk 

(1.9) mes'i n2(iS2 
= NM 

□ 

Lemma 1.3. A monomial of the form f{A, B)a^(3^ appears in deiiW) 
if and only if a monomial of the form g{A, B)a~^[3^ also does. 

Proof: Suppose that f{A, B)a^(3^ is associated to the permutation 



TT. By Lemma p..l| we know that no combination of monomials in the 
expansion cancel. Thus it suffices to display a permutation that pro- 
duces a monomial of the form g{A, B)a~'^l3K We claim that there 
exists a unique permutation vr subject to the following two condi- 
tions: TT picks the entry ((n, m), (n, m)) of iff vr picks the entry 
((n, M — m + 1) , {n, M — m + 1)) , and vr picks {{k,m + 1), (/, m)) iff vr 
picks ((/, M — m + 1), {k, M — m)). It is elementary to verify that vr is a 
permutation that doesn't pick any zeroes. Furthermore, vr and vr pick 
an equal number of elements from each diagonal block, in particular 
an equal number of /3's. vr picks from block (m + l,m) iff vr picks 
from block (M — m + 1, M — m). Thus tt produces a monomial of 
the form g{A, B)a-^(3^ □ 

Notice that, by Lemma |1.2| , if a monomial is of the form /(A, B)a^j3\ 



then f{A, B) is of degree NM—kN—lM. For another pair of exponents 
{k\l'), suppose NM - kN - IM = NM - k' N - I'M. Then {k - 
k')N = (/' - l)M. Since gcd(iV, M) = 1, this implies N\{l' - I). But 
< 1,1 < N. Thus either / = / , in which case k = k , ot I = N 
and I = 0. But if a permutation picks N /3's, it has to be the identity 
permutation. Thus k = and so k = M. We conclude that except for 
the terms a^[3^ and (3^ which have constant coefficients, the degree 
of f{A,B) determines {k,l). 



The following Corollary follows from Lemma's 1.2 and 1.3. 



Corollary 1.1. A monomial with a coefficient of degree d cannot ap- 
pear in det(Vr) unless d is among the following list of numbers: 
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2N, 2N -M, 2A^-2M, ... 

(1.10) NM, NM - M, NM - 2M, NM - 3M, M, 

2NM - 2N, 2NM -2N- M, 2NM -2N - 2M, ... 

2NM - N, 2NM -N-M,... 

2NM 

For 1 < k < M+1, row k of this list contains the numbers {k—l)N—iM 
forO <i< L'^^iJ^J which are all nonnegative. Row M + l + k contains 
(M + k)N — iM for < i < [^^^^^J, i-e. it has the same number of 
entries as row M + 1 — k. 



Notice that the numbers in ( |1.10| ) are situated symmetrically across 



the middle row. We make a definition: 

Definition 1.2. Define i{k) to be the number symmetric to k across 
the middle row in the list ( |1.1(]| ). 

Notice that i{k) = k mod 2, since by definition i{k) — k is 21N where 
/ is the distance between k and the middle row. 

( |1.1CI| ) shows the degrees of coefficients of the monomials. The cor- 
responding 0^13^ are: 

(1.11) a°, oPj3\ 

^-Af+2^^-A./+2^1^^-M+2^2^___ 

It turns out that each of these terms appear in det(Vr) for generic A, B. 
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We inspect the curve det(VF) = when a or P tend to oo. 
and appear in det(PF) with nonzero constant coefficients. has 
coefficient Iln,mB{n,m); we assume that B{n,m) ^ for any {n,m) 
so that this coefficient is not zero. If a is finite and nonzero, it is 
impossible for /3 — > oo on the curve, since the term dominates the 
others. (3 — oo implies a = oo or a = 0, so there are two points of the 
curve at /9 = oo . Call these points P and Q respectively. Similarly, if 
/3 is finite it is impossible to have a = or a = oo. 

The next thing we wish to show is that the behaviour of det(iy) = 
at P and Q is locally identical to the behaviour of + (3^ — and 
Q,-M _|_ _ g these points respectively. More precisely, in the 
normalization X of det(Vr) = 0, there is only one point above either 
of P and Q (these will be denoted by P and Q). This will imply 
div(a) = N{P - Q), and diVoo(/3) = M{P + Q) where a,(5 on X 
mean composition of a, (5 on det(VF) = and projection from X to 
det(l^) = 0. 

To prove this claim, we present algorithms to blow the curve up at 
the points P and Q. These two algorithms are almost identical, so 
we explain the procedure for Q only. Set (5 — ^. So Q is the point 

a — (5 — To make all exponents nonnegative, multiply the curve 
equation by (5^ . The terms in the curve equation are: 



„2M qN 



^2^-2/3^, a2M-2^^-i^ 



, . . . 



[1.12) a^^S^-^, a^/^^'^, a"" p\\a 

a p ,a p , ... 



M 



N 



This curve is singular at Q iff both M and are greater than 1. If 
there is no singularity, there is nothing to prove. So suppose that the 

curve is singular at Q. We blow the curve up at Q, set w = f or w = - 
depending on whether N > M or M > N in the respective order. 
Without loss of generality, we assume that N > M. So we ehminate 
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a from the equation. In the proper transform of the curve, the terms 
are: 

^2M-lpN+M~l ^2M-1 pN+M-2 

^2M-2pN+M-2 ^2M-2pN+M-3 ^2M-2^N+M-A 



:i.l3) ^7^/3^ 



N-M 



Every term except the two boxed terms contain both v and /3 with 
positive exponents (this will be proven below). Hence whether or not 
this blow-up is singular is completely determined by whether or not 
— M or M is 1. Since we assumed that M > 1, it is determined by 
N — M only. If — M is not 1, we blow-up again by setting w = | or 

w = ^ depending on N — M>M or N — M< M, and proceed like 
this until one of the two exponents in the boxed terms is one. 

The following Lemma explains why we can go on, by establishing 
that the two boxes are the deciding terms in every step of the blowing 
up. 

Lemma 1.4. At any stage of the blow-up algorithm, all of the terms 
except the boxed terms (the term on the last row and the term on the 
rightmost column) contain both of the variables of that stage with posi- 
tive exponents. Therefore whether or not there is a singularity is com- 
pletely determined by the exponents of the boxed terms. 

Proof: The proof uses a degree argument. Since we are using the local 
coordinate /3 rather than /3, we use a degree definition slightly different 
from d. Notice, by table ( p..lUD and ( |1.12| ), if we set d{a) = N and 



d{(3) = M, then the boxed terms have degree NM, and all the other 
terms have higher degrees. Whenever a new variable is introduced at 
some step of the algorithm, define its degree naturally as the difference 
of the degrees of the variables that it is a quotient of. Notice that 
all of these degrees are positive because the boxed terms have equal 
total degree, and the comparison of exponents ensures that the new 
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variable is the quotient of the higher degree variable by the lower degree 
variable. Furthermore, notice that, with each proper transform, we are 
decreasing the degrees of all terms in the list by a fixed number. So the 
boxed terms always preserve their significance of being the only terms 
of lowest degree. Now suppose that the variables at some step are v 
and w, boxed terms are v'^ and (so we should have kd{v) — ld{w)), 
and k > I. The new variable of the subsequent step is y = ^ and we 
eliminate w from the equation. Notice that w's are replaced by y's 
with the same exponent, so by induction they survive with a positive 
exponent away from the boxed terms. Contrary to the claim, suppose 
v^w^ is replaced by v^y^ where s < 0. But then i + j — l = s<0. Thus 
i + j < L Since d{v) < d{w) this implies that d{v^w^) < d{w^). But 
this is a contradiction since has minimal degree among the terms. 
□ 

This lemma shows that we may keep blowing up until one of the ex- 
ponents of the boxed terms becomes 1. But then, suppose the sequence 
of variables gotten in the blowing ups is «, f i, "^fc- /5 = implies 

= 0, since in the curve equation at the relevant stage, all but one 
term contains (5. Repeating this, we have ^2 = 0, = 0, = 0. 

But the last one is a nonsingular point, and Vk has a simple zero there. 
Thus there is one point Q over Q in the normalization. Tracing back 
from the last step, from the defining relations of each new variable, it 
is elementary to see that a has a zero of order and (3 has a zero of 
order M at Q. As we remarked, by a similar analysis at P, there is 
one point P above P in the normalization of the curve, a has a pole of 
order iV, and (3 has a pole of order M at P. 

Now, choose a local parameter z at P. By the discussion above 
(a, (3) = {z~^ + + ...), where the Laurent expansions are written 

starting from the lowest degree terms in z . Similarly, at Q, again in a 
local parameter z, (a, (3) — {z^ + + ...). Here, again, the lowest 
degree terms in the expansion are shown. 

We now know the behaviour of the curve at /? = oo. Notice that, 
when P — 0, the curve is defined by the vanishing of a polynomial of 
degree 2M in a. Thus there are 2M points (counted with multiplicity) 
at /9 = 0. Call these points Ri, Si, i = 1, M. We will describe why 
they are indexed like this later. It turns out that the 2M coefficients 
of the equation are algebraically independent, hence almost any set of 
2M points can be gotten by an appropriate choice of A,B. 

In order to calculate the genus, we project X to the /5 axis, and 
use the Riemann-Hurwitz formula. The ramification index at /9 = oo 
is M—1 at each of the two points, so the total is 2(M— 1). Let /(a, (3) ~ 
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■^{a'^^ det{W)). Then the ramification points for finite P values are 
the preimages of points of intersection of / = with det(W^) = 0. We 
need to calculate the degree of this divisor only. Since det(iy) = is 
projective, this number is certainly equal to the intersection number of 
det(iy) = and / = oo. / can be oo only if a or /3 is oo. The only 
points of det(H^) = at this locus are P and Q. Using the given local 
parameters, we inspect the behaviour of / at these points. First of all, 
the terms that appear in / are the a derivatives of the product of a*^ 
and the terms in figure (|1.11|) , i.e. 



(1.14) 

(2M - 1)0=^*^-2, (2M - l)a^^"^f3\ ... 

(2M - 2)a^''-\ (2M - 2)a^''-^(3\ (2M - 2)«2M-3^2^ 

Ma^-\ Ma^-'(3\ Ma^'-^(3^ Ma^'-^(3\ Ma^^-i/5^-\ Ma^'-^(3 

2a^2a^/3^2aV^... 





Therefore at P, the only dominating terms are and 
and in the local parameter z both have Laurent expansions starting 
with These two poles cannot cancel each other, because 

the and [3^ terms cancel in det(Vr), but different constants 
drop in front with a differentiation. 

Similarly, at Q, the term (3^ a^'^~^ dominates by itself, using figure 
1.141 ) again. Its local expansion is ;2-JVAf+{A/-i)Af _^ = ^-iv _^ 



Therefore, the intersection number of / = oo and det = 0, which 
is the total ramification index of X at finite points, equals 2MN — N + 
N = 2MN. 

Now we want to show that the projection X {(3 axis) for a generic 
choice of A,B does not have any other branch points. It is enough 
to show that generic X is nonsingular at finite points. We will prove 
that det(iy) = is nonsingular at finite points for generic A,B. By 
upper semicontinuity theorem, it is sufficient to exhibit one such curve, 
and with no new singularities at oo. Let 77 be a primitive MA^'th root 
of unity. Since gcd(M, iV) = 1, (1,1) generates the additive group 
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Tj/NTj X Z/MZ. Therefore, each (a, a) corresponds to a unique repre- 
sentative (a, a) on the toroidal grid for < a < MN . Let 



a 



(1.15) 5((a,«))=^ 

y4(n, m) = 

Then, first of all 'n.n,m.B{n,m) = rj^ for some s. Any coefficient in the 
curve equation containing an A vanishes. Furthermore, suppose that 
the product C(n,ji),...Xik,jk) — ii)...B{ik, jk) appears as a coefficient 
in some monomial. Since W has toroidal symmetry, C(j^+ijj+i)^...^(j^+ij^+i) 
also has to appear as part of the coefficient associated to the same o^fiK 
Unless f^{ii,ji),...,(ikjk) is the product of all -B's, these two coefficients are 
different. When i?'s assume the values in (|1.15| ), the ratio of the second 
coefficient to the first is rj ^ But Y.n=i V^"" = for any k < NM. Thus 
the entire sum vanishes for such k. We deduce that the curve equation 
for this assignment of values is: 



(1.16) «2M ^ ^A/^TV ^ _ Q 

It easy to see that (|1.16|) has no singular points other than P and Q 
with the multiplicities described, so the claim is established. 

Therefore, for the generic curve, the total branching is the calculated 
minimum, 2NM + 2M — 2, not more. By Riemann-Hurwitz formula. 



2NM + 2M -2 

, 9 = 2M + 1 

(1.17) ^2 

= {N -\)M 

Next, we look at the solutions \E' of = 0. By linear algebra, the 
components "^{kj) of should obey 

,^ ^(A;,/) _ det(iy(^'^)'('='')) 

1 1.18 ) 



^(a,6) det(W^(^'J)'(»''')) 

for any row (recall that we use a pair of indices to indicate a 

row or column). Here, denotes the minor of W obtained by 

removing row and column {k,l). 

Once and for all, we normalize \E' such that \1'(0, 0) = 1. Then \E'(A;, /) 
is a meromorphic function on X — {P, Q} hj (|1.18|) . We first want to 
compare the behaviour of the "^{k, /)'s at the /5— infinite points P and 
Q of X. We look at ^|^^7j^, and ^^^^jy-- We make a careful choice of 
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(i, j) for these analyses: Use (z,j) = {k,l + 1) for both comparisons. 
Repeating the formula above, we have: 



det(Vr('=''+i)'(*=+i'')) 
det(iy('=''+i)'(^'^)) 



Although it is definitely difficult to find the determinants on the right 
hand side explicitly, it is much easier to determine their behaviour at P 
and Q. Notice that the element {{k,l + 1), {k, /)) of W is A{k, I). This 
immediately implies that the expansion of det(iy('=''+i)'(^'')) contains an 
a°-(3'' term iff the total coefficient of q;"/5^ in det(Vr) contains A{k,l). 
The coefficients in det(iy) are toroidally symmetric, therefore we may 
deduce that this happens iff the coefficient of contains any A. 
From this observation, one can easily make a list of the a"/3^ that 
appear in det(W^*^'^''"'"^'''*-'^''^), but the following result will be enough for 
us: 

Lemma 1.5. In the expansion o/det(iy*-'^''"'"^-*'*''^''''), the terms , 
and (3^ are absent, the terms with coefficients of degree 1 and degree 
i{l) (see definition \1.2J are present. 

Proof: The absence result is clear. The coefficients of a^^, or 
13^ never contain an A. For the presence result, note that 1 and i{l) 
are both odd numbers. Therefore a monomial consisting in only 5's, 
which evidently is of even degree, cannot give them. □ 

We make a similar analysis for det(iy*^'^''"'"^'''''''''"'"^^). The element 
((/c, / + 1), (/c, / + 1)) of FT is a diagonal entry. Hence the terms that 
appear in this determinant are precisely those whose coefficients in 
the expansion of det(iy) contain a diagonal entry, and consequently 
contain a [3. Unless the ((/c, / + 1), (fc, / + 1)) term is itself [3 (i.e. unless 
/ = 0) the complete list of terms in the expansion of (iei{W'^^'^^^'>''^'^'^^^^) 
is: 



:i.i9) 



and 



:i.2o) 



^(fc+i,/) 

v|>(fc,0 



^(A:,/ + l^ 
v&(A;,/) 
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(1.21) a'p\ a'P^ 

The last determinant that we need to look at is det(W^*^'^''+^'''*^''"'"^'')). 
Similarly, look at terms whose coefficients in det(iy) contain {{k,l + 
1), {k + 1,/)). This is an off-diagonal 1. Again, we are going to state 
what we need: 

Lemma 1.6. In the expansion of det{W^^'''~^^^'^''^^'^^) , the a'^^ and de- 
gree terms are absent, the and degree i{2) terms are present. 

Proof: The absence of the o;"*^ term and the presence of term 
are clear. By the symmetry of Lemma , an off-diagonal 1 appears 



in the degree z(l) coefficient iff a appears in the degree 1 coefficient. 
This is impossible since d{B) = 2. Similarly, a 1 appears in the i{2) 
coefficient iff a B appears in the degree 2 coefficient. This always 
happens: Pick B{r,s), and now set all other A,B = Om W. Now 
there are two nonzero elements in every row and column except the 
ones that B{r, s) belongs to. The only nonzero element left in column 
(r, s + 1) is ((r — l,s + 2), (r, s + 1)). Pick this. Continuing, pick 
all {{r — t, s + t + 1) , {r — t + 1, s + t)) until this sequence comes to 
((r — 2, s + 1), (r — 1, s)). This point appears before {{r, s + 1) , {r + 1, s)) 
in the sequence because gcd(A^, M) = 1. Pick all the diagonal elements 
for the remaining rows and columns. This gives a valid permutation of 
the type we want. □ 

By using Lemmas [L^, and figure (|1.21| ), we can compare the 
behaviour of different Z) at P and Q 

Proposition 1.1. has one pole at P and one zero at Q. 

Proof: At P, the dominating term of det(W^(''''+^)'(*''+^'')) is a*^. Ex- 
panded in the local parameter, this is z~^'^^ + .... The dominating 
term of det(Vr('=''+i)'('='')) is the de gree 1 term, and the local series is 
^-MN+i _|_ ____ Thus the claim about P is established. At Q, the domi- 
nating terms have degree i{2) and degree coefficients respectively. 
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Thus one gets z •^•^+2 -|- ... and z ^^+1 + ... in the local parameter 
respectively. This finishes the proof of the claim about Q. □ 

Proposition 1.2. ^^'^^^ has one pole at P and one pole at Q. 

Proof: At P, the dominating terms of dct(l^(''''+^)'(''''+^)) and det(W^('''^+^) 
are (3^ and the term with degree 1 coefficient respectively. Local ex- 
pansions give z"^^ + ... and z~^^^'^ + .... And at Q, the dominating 
terms are and the term with degree coefficient respectively. 
These give z~^^ + ... and z~^^^^ + ... respectively. □ 

So we have obtained that ratios of ^' obey this fixed structure of 
poles at /? = 00. Now we look at /3 = 0. Notice that when /3 = 0, 
det(Vr) splits as det(A:(l))... det(X(M)). Likewise, when /3 = the 
^Q{(W{k,i+i)m)) and det(Vr('=''+^)'(*^+i'')) split as a product of block 
determinants. Observe that in these block expansions, all det(X(m)) 
appear except for det(X(/)). The remaining A^ — 1 x A^— 1 determinant 
has nonvanishing determinant at points Si of the curve for generic 
A,B. Therefore, these two determinants have the same zeros on the 
curve at /9 = 0. This shows that ^(A;, /) and \E'(fc + 1, /) have the same 
divisor at /3 = 0. On the other hand, we saw that det(L('^'^+^)'*^''''"'"^)) 
vanishes identically at /? = 0. Therefore \E'(/c, / + 1) has two more zeros 
than '^{k,l) at /3 = 0, namely the two solutions for a of det(X(/)). 
Call these two points Ri, Si. 

The line bundle mentioned in the correspondence comes from a di- 
visor V encoding the totality of remaining poles of \l/'s: 

Definition 1.3. Let V be a minimal effective divisor such that 

m 

V + div{-^{n, m)) + (n + m)P + (m - n)Q - + 5'^) > 



i=l 



for all {n, m). 



We remark that such a finite V exists: Suppose the condition is 
satisfied for < n < A^ and < m < M, then it is automatically 
satisfied for all positive (n,m). This happens since ^(n + kN,m + 
IM) — a''(3'''i/{n,m). The new divisors introduced by the a's and /9's 
are exactly taken care of by the extra points in definition of T>. Define 
the fine bundle C — 0{V). Furthermore, set 



[1.22) Ln,m = C{{m + n)P + (m - n)Q - J^i^i + Si)) 



i=l 



The choice of T> says ^(n, m) is a section of Cn,m- 
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Now, fix m,n, and select k large enough so that deg{T>) + k > 2g — 2. 
For this particular k, Riemann-Roch theorem gives 

(1.23) h'^iCnAkP)) = deg{Cn,m) + k-g + l 
because h}{Cn,m{kP)) = by the choice of k. 

Proposition 1.3. ^(n + k,m) is a section of Cn,m{kP) but is not a 
section of Cn,miik — 1)P) 

Proof: \I'(n + k,m) is a section of Cn,m{kP) by definition. Suppose 
it is a section of Cn,m{{k — 1)P). Then by the pole comparisons at 
P, we can replace k by any nonnegative number, and the hypothesis 
remains true. But then \E'(?T,,m) is a section of Cn,mi—P), and by pole 
comparisons again this holds for all [n, m). In particular ^1^(0, 0) = 1 is 
a section of C{—P), therefore V — P is effective. But these imply that 
we can replace V hj V — P, which contradicts the minimality of V. □ 
From the proposition we deduce that 

(1.24) h%C^,rn{kP)) = deg{V) + k - g + 1 

for all k, in particular h^{Cn,m) = degiV) — g + 1. Indeed, increasing k 
by one always gives a new section, therefore, should be zero at each 
step, because it is eventually zero. 

The divisor V is of degree g. For a proof of this, we refer the reader 
tofCg. 

This concludes the construction of the algebro-geometric data from 
the difference operator. We now summarize the construction of the KP 
data from the algebro-geometric data. 

Assume A" is a smooth curve of genus g^ P, Q, Ri, Si points on X 
such that div(a) = N{P-Q) and div(/3) = MiP + Q^-'^fi^iRi + Si), 
L a line bundle of degree g ot^ X such that 

m 

(1.25) h^{X, £((n + m - \)P + (m - n)Q - ^(i?, + 5,))) = 
for all m, n. Let 

ra 

(1.26) Ln^m = C{{n + m)P + (m - n)Q - ^(i?, + Si)) 

i=l 

as before. Then, 
Proposition 1.4. /i°(A',£„^m) = 1 for eachn,m. 
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Proof: By Riemann-Roch theorem, h^{X,Cn,m) > 9 + ^ ~ 9 = 1- 
On the other hand, it is impossible to have strict inequahty because 
of the regularity assumption. (Removing one P prohibits at most one 
section.) □ 

Say ^(r2,m) is a nonzero section of Cn,m- 

Proposition 1.5. The set of sections {\E'(?2 + 1, m), \E'(n, m), \E'(?2 — 
l,m), \E'(?2, m + 1)} is linearly dependent. 

Proof: Notice that all four of them are sections of Cn,m{P+Q)- On the 
other hand, by Riemann-Roch theorem, we get h^{Cn,m{,P + Q)) > 3, 
and by the regularity assumption again, equahty has to be realized. □ 

The ratio fn,m = '^^ll^^-j has a pole of order precisely n + m at P. 
Suppose 2 is a uniformizing parameter at P. We can normalize \E'(n, m) 
so that fn,mZ^^'^ = 1 at P. We can also normalize the meromorphic 
functions a, (3 such that a = + P = z~^^ + .... 

After these normalizations it is elementary to show that the linear 



dependence between the four sections in Prop 1.5 has to be of the form 



(1.27) 

\l/(n + 1, m) — \E'(n, m + 1) — A{n, m)'^{n, m) — B{n, m)\l'(n, m) = 

for some functions A, B on Z/NZ x Z/MZ. 

Thus we get back to A,B, having started from the algebro-geometric 
data. Next, we wish to calculate the flow equations explicitly by com- 
puting the effect of moving C in the Jacobian in certain directions. 



2. The flow equations 

2.1. Some functions on Z/NZ x Z/MZ. As before, suppose that 
N,M G Z, and gcd(A^, M) = 1. Let S denote the set of functions 
/ : Z/NZ X Z/MZ ^ {-1, 0, 1}. 

Proposition 2.1. There is a unique function n in S that satisfies the 
following conditions: 



(i)K(0,0) - k(1,-1) = -1 
(ii)K(0,l) - /«(1,0) = 1 
^ ' (m)K(-l,0) - k(0,-1) = 1 

(if)K(-l, 1) -k(0,0) = -1 

and except for these four values of {i,j), n{i — 1, j + 1) = 
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Proof: Uniqueness is easy to prove, because if two such functions 
exist, their difference has to be a constant. But by (i) and (iv) the only 
possibihty for k(0, 0) is 0. Therefore the constant is zero. 

To prove existence, we note that since gcd(A^, M) = 1, (—1, 1) is a 
generator. Look at the sequence 



(2.2) (-l,l),(-2,2),...,(-a,a),... 

We will distinguish the two cases below: 

1) Suppose in the sequence (|2.2|), (1, 0) appears before (—1, 0). Then 
we declare 



«:(-!, l) = «:(-2, 2) = ... = «:(1,0) = -1 
^ ' ' -1) = k{2, -2) = ... = 0) = 1 

and K,{a, 6) = if (a, b) is not in these lists. 

2) Suppose (—1,0) appears before (1,0). We declare 

1) = 2) = ... = ^(0, -1) = -1 

^ ■ ' k(1, -1) = k{2, -2) = ... = k(0, 1) = 1 

and K{a, 5) = if (a, b) is not in these lists. 

One can check that k satisfies the conditions that we asked for. □ 
Remark: One may wonder which case happens when. It turns out 
that the deciding quantity is the parity of the number of steps in the 
Euclidean algorithm for the ordered pair {N, M). In particular we have 
alternate cases for {N,M) and {M,N). 

Note that K{n,m) = —K{—n, —m). 

We introduce another function p in S: 

Definition 2.1. 

(2.5) p(n, m) = K{n + 1, m) + /t(n, m) + 5(^n,m),{o,o) - 5(n,m), (-1,0) 

(2.6) (j){n,m) = —p{—n — 1, — m) — p{~n, —m) 

{5 is the Kronecker delta function, i.e. 5x.y = 1 if X = F, and 
otherwise, etc.) 

Then the following holds: 

p{n — 1, m + 1) — p{n, m) ={K,{n, m + 1) — K,{n + 1, m)) 

+ {K,{n — 1, m + 1) — K{n, m)) 

+ f (n-l,m+l),(0,0) — f (n_i,m+l),(-l,0) 
— 5(n,m), (0,0) + ^(n,m), (-1,0) 
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Therefore p is the unique function in S satisfying the following condi- 
tions: 

(,)p(_2,0)-p(-l,-l) = l 

(n)p(0,l)-p(l,0) = l 
^ ■ ^ (m)p(-l,0)-p(0,-l) = -l 

(i^;K-l,l)-p(0,0) = -l 
and p{n — 1, m + 1) = m) for all other (n, m). 

2.2. Derivation of the discrete KP equations. We would hke to 

compute what happens to A and B when we keep X and the points 
fixed, but deform the line bundle C linearly in the Jacobian. By choos- 
ing the direction of deformation carefully, the computation of the first 
few evolution equations for A and B becomes feasible. 

X is mapped inside its Jacobian by the Abel-Jacobi map .We are 
going to deform the line bundle in the tangent direction to the curve 
at P. Assume again that 2; is a local parameter at P. For small t e C, 
define 



(2.9) f~'n,m,t — ^n,m{P ^) 

where T = z~^{t). First, notice that h'^{Cn,m,t{~-P)) = by upper 
semicontinuity for small t, since hP{Cn,m{~P)) — 0- Thus hP{Cn^m,t) — 
1. Suppose ^(n, m, is a nonzero section. 

Proposition 2.2. The set {^{n,m,t),^{n,m),'^{n -\- l,m)} is lin- 
early dependent. 

Proof: This follows from the observation that all three are sections of 
^^n,m,t{T) and [Cn,m,tiT)) = 2 by Riemann-Roch theorem. □ 
Thus suppose 

(2.10) 

^(n, m, t){z) — ao{n, m, t)'i>{n, m){z) + ai{n, m, t)'^{n + 1, m){z) 

where ao and Oi are holomorphic in t. Now, ^(n, m, t) vanishes at z — t 
by definition. Therefore 



(2.11) ao{n, m, t)^(n, m){t) + ai{n, m, t)^'(n + 1, m){t) = 

Now we can expand this expression in powers of t and calculate the first 
few terms of the power series expansions of the holomorphic functions 
tti in terms of the Laurent coefficients of the function '^^^^^^ whose pole 
behaviour at P we know. Let's give a name to these coefficients. 
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First define 



(2.12) f{n,m) , , 

^ ' M , ; ^(0,0) 

/(n, m) has a pole of order n + m at P. Tlien let diin^ m) be sucli tliat 



(2.13) f{n, m) = z"""™ + di{n, 771)^-""™+^ + d2(n, m)2-"-'"+2 + ... 

Tlien since 
(2.14) 

\I/(n + 1, m) = \l/(n, m + 1) + A(?2, m)\l/(n, m) + B{n, m)\l/(n — 1, m) 
we liave 

(2.15) di{n,m) — di{n,m + 1) = A{n,m) 

Now we can expand ( |2.11|) in powers of t. We define one more piece of 
notation for tliis purpose: say ai(n, m, t) = aij{n, m)tK 

One sees immediately that ai i(?t,, m) = —1, therefore we can normal- 
ize ^(n, m, t) in a neighborhood of P, replacing it by \E'(r;,, m, t)t/ai{n, m, t). 

Then if we compute the coefficient of t~^~'^~^^ we get 

(2.16) 00,1(^^5 nT-) = di{n + 1, m) — di{n, m) 
On the other hand, one can expand the relation 
(2.17) 

\E'(n + 1, m, t) = \E'(n, m + l,t) + A{n, m, t)\I^(n, m, t) + i?(n, m, t)\&(n - 1, m, t) 

using equation (|2.1CI|) . One can further write all \E'(fc,/) in terms of 
\E'(n + i,m) by reducing I to m using the difference relation repeatedly. 
Finally, one gets a time dependent linear relation between \E'(fc, m), 
where m is fixed. But these all have different order poles at P as 
we have seen, so they are linearly independent. Therefore many time 
dependent quantities (that one may calculate from the last equation) 
ought to vanish. The coefficients of t in front of \E'(n, m) and \E'(?7,— 1, m) 
give evolution equations for A and B: 

(2.18) 

A{n, m) = ai,i(n, m + l)B{n + 1, m) + (ao,i(n, m + 1) — ao,i('^) m))A{n, m) 

— 0-1,1(71 — 1, m)B{n, m) 
B{n, m) = (ao,i(n, m + 1) — ao,i(n — 1, m))B{n, m) 
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Writing a's in terms of ci's, equations become 
(2.19) 

A{n, m) = B{n, m) — B{n + 1, m) + {di{n + 1, m + 1) — di{n, m + 1) 

— di{n + 1, m) + di{n, ■m))A{n, m) 
B{n, m) = {di{n + l,m + 1) — di{n,m + 1) — di{n, m) 

+ di{n — l,m))B{n, m) 



rfi(0, 0) = by definition. Now, in tlie light of ( |2.15| ), the terms in both 
parentheses are hnear polynomials in ^'s. 

(2.20) 

A{n, m) =B{n, m) - B{n + 1, m) + (^ g(^n,rn),{k,r)A{k, l))A{n, m) 

B{n, m) =(^ h(^n,m)XKi)A{k, l))B{n, m) 
By toroidal symmetry, 
(2.21) 

for any (a, 6) G Z/A^Z x Z/AfZ. Same result holds for h. Thus, with 
some abuse of notation, we may write 

(2.22) 

A{n^ m) =B{n, m) — B{n + 1, m) + (^^ g{k — n,l — m)A{k, l))A{n, m) 
B{n, m) =(^^ h{k — n,l — m)A{k, l))B{n, m) 

We subtract the terms of the first equation of ( |2.22| ) in parentheses for 
(n, m) from the same terms for (n + 1, m — 1): 

(2.23) 

{J29ik-n-l,l-m + l)A{k,l)) - {J^gik - n,l - m)A{k,l)) = 
{di{n + 2, m) — di{n + 1, m) — di{n + 2,m — 1) + di{n + l,m — 1)) 

— {di{n + 1, m + 1) — di{n, m + 1) — di{n + 1, m) + di{n, m)) 
= A{n + 1, m) — A{n, m) — A{n + 1, m — 1) + A{n, m — 1) 

And we do the same for the second equation: 
(2.24) 

(^h{k - n - 1,1 - m + l)A{k, h{k - n,l - m)A{k, I)) = 

{di{n + 2, m) — di{n + l,m) — di{n + l,m ~ 1) + di{n, m — 1)) 

— {di{n + 1, m + 1) — di{n, m + 1) — di{n, m) + di{n — 1, m)) 
= A{n + 1, m) — A{n, m) — A{n, m — 1) + A{n — 1, m — 1) 
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Equating the coefficients for all {k, I) we obtain 

(i)^(0,l)-^(l,0) = l 

{n)9{-lA) - 9M ^ -I 
^ ■ ^ (m)^(0,0)-^(l,-l) = -l 

{iv)g{-l,0)-g{0,-l) = l 

and g{n + 1, m — 1) = g{n, m) for any other (n, m). 
and 

(i)/i(0,l) - /i(l,0) = 1 

- /?.(0,0) = -1 

(m)/?.(-hO) - /?.(0. -1) = -1 

(«f)/i(-2,0) - /i(-l,-l) = 1 

and + 1, m — 1) = m) for any other (n, m) 

These are precisely the relations for k and p, so we have proven: 

Proposition 2.3. T/ie equations of evolution are 
(2.27) 

74(n, m) =B{n, m) — B{n + 1, m) + — n,l — m)A{k, l))A{n, 

B{n, m) —C^^ p{k — n,l — 'm)A{k, l))B{n, m) 



(2.26) 
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